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ABSTRACT 



Recent XMM-Newton observations have measured the hard (2.5-12 keV) X-ray 
spectrum of the well-known Seyfert galaxy NGC 4151 with a signal-to-noise unprece- 
dented for this source. We find that a spectral model, developed to fit previous Bep- 
poSAX and ASCA observations of NGC 4151, provides an excellent description of 
the XMM-Newton EPIC data. The results support the view that it is the level of 
the continuum that is the main driver of the complex spectral variability exhibited 
by NGC 4151. We focus on the iron K features in the NGC 4151 spectrum. There 
is no requirement for a relativistically broadened iron Ka line, in contrast to several 
earlier studies. The iron Ka line profile is well modelled by a narrow Gaussian, the 
intensity of which varies by ~25% on timescales of about a year. There is also a strong 
suggestion that the cold media present in the active nucleus of NGC 4151 have an iron 
abundance that is at least twice the solar value. 

Key words: galaxies: active - galaxies: Seyfert - X-rays: galaxies - galaxies: NGC 
4151. Galaxies 



1 INTRODUCTION 

The Seyfert 1 galaxy NGC 4151 harbours one of the bright- 
, est Active Galactic Nuclei (AGN) accessible in the X-ray 
band and has been extensively studied by all major X-ray 
missions. This observational focus has revealed the X-ray 
spectrum of NGC 4151 to be a complex mixture of emis- 
sion and absorption components, originating at a variety 
of locations from the innermost parts of the putative ac- 
cretion disk, out to the extended narrow-line region of the 
galaxy. The X-ray spectrum emanating from the active nu- 
cleus is dominated by an intrinsic X-ray to 7-ray contin- 
uum probably produced by the thermal Comptonization of 
soft seed photons (e.g. Zdziarski et al. 1996, 2000). There 
is also evidence in the X-ray spectrum for continuum repro- 
cessing features, in the form of a Compton reflection con- 
tinuum and emission/absorption features imprinted by iron 
(Yaqoob et al. 1993; Zdziarski et al. 2002; Piro et al. 2002). 
Below ~ 5 keV the hard continuum is strongly cut-off by 
photoelectric absorption in a substantial (A?jj~10 23 cm -2 ) 
line-of-sight gas column density (Holt et al. 1990; Yaqoob 
et al. 1993; Weaver et al. 1994b; Schurch & Warwick 2002). 
Additional soft X-ray emission becomes apparent below ~2 
keV and dominates the spectrum below ~1 keV (Weaver et 
al. 1994a,b, Warwick, Smith & Done 1995). This soft emis- 
sion has been spatially resolved by Chandra(Og\e et al. 2000; 



Yang, Wilson & Ferruit 2001) and is co-spatial with the O 
III ionization cones identified in Hubble Space Telescope ob- 
servations (Ogle et al. 2000) implying a strong association 
between the soft X-ray emission and the extended narrow- 
line region (ENLR) of is galaxy. 

Some earlier studies, based on ASCA observations, con- 
cluded that the profile of the iron Ka emission line is com- 
plex and may be composed of an intrinsically narrow com- 
ponent plus a relativistically broadened line feature (Yaqoob 
et al. 1995; Wang et al. 2001). Furthermore, the broad line 
profile has been reported to be variable on timescales of 10 4 s, 
corresponding to an emitting region of <0.02 AU, suggesting 
an origin close to the supermassive black hole presumably 
in the inner regions of a putative accretion disk (Wang et 
al. 2001). However, there are alternative interpretations of 
the broadband ASCA spectra that do not require the pres- 
ence of an extremely broad iron line in NGC 4151 (Schurch 
& Warwick 2002; Takahashi et al. 2002). 

In this paper we focus on the iron K features present 
in the extremely high signal-to-noise spectra of NGC 4151 
recorded by the EPIC CCD cameras on XMM-Newton. Our 
analysis makes use of the spectral "template" model, de- 
veloped by Schurch & Warwick (2002) to interpret earlier 
ASCA and BeppoSAX observations. 
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2 THE XMM-Newton OBSERVATIONS 

NGC 4151 was observed by XMM-Newton three times dur- 
ing the period December, 21 si -23 rd 2000 (orbit 190). In 
all three observations the EPIC MOS and PN cameras 
(Turner et al. 2001; Striider et al. 2001) were operated with 
the medium filter in place. The first observation (~33 ks) 
was carried out with both CCD systems in Small Window 
Mode (SWM), whereas the remaining two observations (~23 
ks and ~63 ks respectively) employed Full Window Mode 
(FWM). The events recorded in FWM were screened with 
the XMM-Newton Science Analysis Software (SAS v5.2) to 
remove known hot pixels and other data flagged as bad. 
The SWM observation was processed with a local develop- 
ers version of the SAS, which included a patch with the most 
up-to-date correction for the PN charge transfer inefficiency 
(this correction is now included in the latest public release 
version of the SAS, v5.3). The data were processed using the 
latest CCD gain values and only X-ray events corresponding 
to patterns 0-12 in the MOS cameras and 0-4 in the PN cam- 
era were accepted. An examination of the temporal variation 
in count rate from the full available field (but excluding the 
bulk of the NGC 4151 contribution) in all three observa- 
tions revealed a single background flaring event during the 
first observation, which was screened from subsequent anal- 
ysis. The resulting total effective exposure times were 110 
ks and 91 ks for the MOS and PN instruments respectively. 

An investigation into the impact of pile-up on the ob- 
servation showed that the effect was negligible in almost all 
cases, largely due to the low flux state of NGC 4151 during 
the observation period. The exception to this is the FWM 
MOS observation, which shows marginal pile-up above 8 
keV. The best-fit spectral parameters and \ 2 values quoted 
in this paper are for spectral fits in which the MOS FWM 
data above 8 keV have been excluded. 



3 THE X-RAY LIGHT-CURVE 

A source light-curve was extracted for the MOS cameras 
from each individual, screened observation. Source counts 
were taken from within a 100"xl00" box for the SWM ob- 
servation and a 2' radius circle for the FWM observations. 
Background subtraction was not applied since the back- 
ground contributes no more than a few percent of the total 
signal. 

The 2-10 keV light-curve (Fig. 1, Upper Panel) shows 
slow drifts in the X-ray flux of NGC 4151 on timescales 
longer than ~10 4 s, behaviour which is very characteris- 
tic of this source (e.g. Yaqoob & Warwick 1991; Yaqoob et 
al. 1993). The largest flux variation seen was an increase 
of ~25% over the last ~2xl0 4 s of the observation. Previ- 
ous observations of NGC 4151 have revealed complex spec- 
tral changes associated with long-timescale flux variations 
(Schurch & Warwick 2002). However, the 2-5 keV/5-10 keV 
softness ratio measured during the XMM-Newton observa- 
tions (Fig. 1, Lower Panel) remains fairly constant, with just 
a hint of spectral softening towards the end of the observa- 
tion correlated with the flux increase noted earlier. 







i ■ i 1 i 


i i i 


— 1 1 1— 



2x10 4 4x10 4 6x10 4 8x10 4 10 5 1.2x10 5 



Time (s) 

Figure 1. Upper Panel: The 2-10 keV light-curve from the co- 
added XMM-Newton EPIC MOS detectors. The light-curve is 
binned to a resolution of 500 seconds and covers the three in- 
dividual observations (separated by the data gaps) . Lower Panel: 
The variation of the 2-5 keV/5-10 keV softness ratio during the 
observations. 

4 SPECTRAL ANALYSIS 

In the following analysis we consider the background- 
subtracted spectrum of NGC 4151 averaged over the three 
XMM-Newton observations. For the MOS spectra, source 
counts were extracted from the same source regions as em- 
ployed for the light-curves. In the FWM observations, the 
background was taken from a region offset from the source 
response but on the same central chip (MOS 1 and 2 cam- 
eras). The FWM background was also used to correct the 
SWM observation. For all the MOS observations the total 
background amounted to less than 3% of the source plus 
background count rate. 

For the PN camera we used a 2' radius source cell for 
the FWM observations with an offset background region lo- 
cated on the same chip as the source image. For the SWM 
observation we derived the source spectrum using a 4'x2'box 
and the background from a similar size region offset from the 
source but on the same chip. In this latter case the contam- 
ination of the background region by the wings of the source 
response resulted in a background count rate amounting to 
~8% of the source plus background signal. The loss of source 
flux in the various cases is accounted for both in the spectral 
fitting and in the quoted fluxes. 

For convenience the spectra from the two FWM obser- 
vations were co-added to produce a single source spectrum 
and a single background spectrum for each instrument. The 
resulting spectra were binned to a minimum of 20 counts per 
spectral channel, in order to apply \ 2 minimisation tech- 
niques in the spectral fitting process. Here we restrict the 
following analysis to the 2.5-12 keV bandpass. A study of 
the complex soft X-ray spectrum of NGC 4151 based on 
XMM-Newton EPIC and RGS data is in progress and will 
be presented elsewhere. 

We adopt the spectral "template" model described by 
Schurch & Warwick (2002), which includes the following 
emission components: 

(i) An absorbed power- law continuum with a fixed photon 
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Figure 2. Upper Panel: The 2.5-12 keV X-ray spectra measured 
by the EPIC instrument on XMM-Newton. The MOS data, shown 
here in black (MOS 1) and red (MOS 2), are from the SWM ob- 
servation. The PN data, shown here in green, are from the FWM 
observation. The solid line shows the best-fitting spectral "tem- 
plate" model (Table 1, Model 1). Lower Panel: The ratio of the 
data to the model prediction. 



index, F=1.65, and a high energy break at 100 keV. The 
normalization of the power-law, Ai, is a free parameter in 
the modelling. 

(ii) A neutral Compton reflection component (repre- 
sented by the XSPEC model pexrav; Magdziarz & 
Zdziarski 1995) with only the reflection scaling factor, R, as 
a free parameter. Specifically, the normalisation and spectral 
index of the illuminating component were fixed at the same 
values as those for the underlying power-law continuum. In 
addition cos i was set at 0.5 and the metal abundance in the 
reflector was taken to be the solar value. 

(iii) An iron Ka emission line with a Gaussian profile, 
with the line intensity Ikq, the line energy Ekq and the 
intrinsic line width cjkc all free parameters. 

The absorption of the power-law continuum (which is 
known to be complex in NGC 4151) is represented as the 
product of two absorption components, namely a partially 
photoionized (i.e. warm gas) component with column den- 
sity ^H,marm and a cold gas column Nn,coid- The state of 
the warm absorber is governed by the ionization parame- 
ter £ defined as £=Li 0n /nr 2 , where Li 0n is the source lumi- 
nosity in the 0.0136-13.6 keV bandpass (erg s -1 ), n is the 
number of hydrogen atoms/ions in the gas (cm -3 ) and r 
is the distance from the central source to the inner edge of 
the warm cloud (cm). The value of £ is a free parameter 
in the model, but the column densities of the two absorp- 
tion components are fixed at the values used by Schurch 
& Warwick (2002), namely N HiUjarm =2.4x 10 23 cm -2 and 
^H,coid=3.4x 10 22 cm" 2 . For further details of the pho- 
toionization modelling see Schurch & Warwick (2002) and 
Griffiths et al. (1998). The adopted spectral model also in- 
cludes absorption arising in the line-of-sight column density 
through our own Galaxy (N£r,Goj=2xl0 20 cm -2 ), applied 
to all the emission components. 

In summary, the above spectral model requires a total 
of just six free parameters relating to the continuum and 
reflection normalisations, the iron Ka line properties and 
the ionization state of the warm column. However, in the 



Table 1. The best-fitting spectral parameters 



Parameter 


Model 1 


Model 2 


R 

log £ 
A pl a 

0~Ka b 
l Ka C 
^K,Edge b 
T~K,Edge 


2 634+0 003 

1 r,o + 0.03 

L - Z °-0 01 
g ogo+0.003 
o.oyo_ 002 

03Q + '005 

i o 9 +0.03 
1 - oz -0.03 


1 oq+0.14 

2 624 +o: ° 04 
z - uz -0.006 
i ,1+0.03 
^■ o± —0 02 
6 393+0.003 

n 33+0.006 

1 26 +0 - 04 
i ' ZD -0.04 
7 ii+O.Ol 

'• 11 -o.oi 

12+ - 02 
u - ±z -0.02 


x 2 


5049 


4859 


d.o.f 


4756 


4754 



a 10~ 2 photon koV -1 cm~ 2 s 
b keV 

c 10~ 4 photon cm -2 s -1 



simultaneous fitting of the 5WMand FWM spectra for each 
detector (MOS 1, MOS 2, PN), we also employ five constant 
factors allowing for up to a 5% relative flux scaling between 
the different spectral datasets. Although the resulting best- 
fit was formally unacceptable, with x 2 =5193 for 4756 de- 
grees of freedom (d.o.f), the remaining residuals were largely 
at the few percent level, close to the limit of the calibration 
uncertainty of the XMM-NewtonEPIC cameras. Including a 
1.5% systematic error to allow for these calibration uncer- 
tainties (S. Sembay, private communication) reduces the x 2 
to 5049. By way of illustration, Fig. 2 shows three of the 
six available EPIC spectra datasets along with the overall 
best-fitting model and the data/model fit residuals. The cor- 
responding best-fitting parameter values are listed in Table 1 
as Model 1. The quoted errors (here as elsewhere in this pa- 
per) are for a 90% confidence level as defined by a Ax 2 =2.71 
criterion (i.e., assuming one interesting parameter). 

Remarkably, the spectral template model provides an 
excellent description of the XMM-Newton spectra with only 
minor adjustments to a limited number of free parameters. 
Given the very high quality of the EPIC spectra (the PN and 
the combined MOS 1/2 spectra contain in excess of 340,000 
and 240,000 counts in the hard X-ray band, respectively!), it 
would be surprising if further subtle spectral features were 
not present in the data. In fact, inspection of the data/model 
residuals (see Figure 2, Lower Panel) do suggest a systematic 
deficit of counts between 7 and 8 keV. This is interpreted as 
a need for extra absorption over and above the current model 
prediction, which in terms of the spectral modelling can be 
accommodated by the addition of an extra absorption edge, 
with the edge energy and optical depth as free parameters. 
Also at this stage a contribution from (neutral) iron K/3 
emission is included; for this purpose a Gaussian line was 
added at a rest frame energy of 7.058 keV, with the same 
intrinsic width as the iron Ka line and at 10% of the Ka 
intensity (i.e. the appropriate branching ratio). 

The result of these modifications was a significantly im- 
proved fit (x 2 =4859 for 4754 d.o.f). The model fit to the data 
is shown in Fig. 3 and the best-fitting parameter values for 
this revised model are listed as Model 2 in Table 1. The null- 
hypothesis probability of this latter model (0.14) indicates 
that this is a reasonable representation of the spectral data 
to within the limits of the instrument calibration. 
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Figure 3. Upper Panel: The 2.5-12 keV X-ray spectra measured 
by the EPIC instrument on XMM-Newton. The MOS data, shown 
here in black (MOS 1) and red (MOS 2), are from the SWM 
observation. The PN data, shown here in green, are from the FWM 
observation. The solid line shows the best-fitting spectral model 
(Table 1; Model 2). Lower Panel: The ratio of the data to the 
model prediction. 



5 DISCUSSION 

On the basis of the models in Table 1, the absorption cor- 
rected flux of NGC 4151 measured by XMM-Newton was 
5.8xl0 -11 erg s _1 (2-10 keV). This represents a very faint 
state of the source as judged against previous extensive X- 
ray monitoring by missions such as EX OS AT, Ginga, Bep- 
poSAX and ASCA (Yaqoob et al. 1993; Piro et al. 2002; 
Schurch & Warwick 2002). 

A number of points follow from a detailed comparison 
of the parameter values in Table 1 (Model 2) with those 
quoted in Schurch & Warwick (2002). For example, the ion- 
ization parameter pertaining to the warm absorbing medium 
remains high (log £=2.62) despite the decline in the level of 
the hard power-law continuum by, for example, almost a fac- 
tor 4 compared with the January 1999 BeppoSAX observa- 
tion. This may be further evidence for long recombination 
timescales in the warm absorber as discussed in detail by 
Schurch & Warwick (2002). In actual fact, the value of the 
ionization parameter derived from the XMM-Newton obser- 
vation falls at the high end of the range recorded in the 
ASCA long- look observation of May, 2000 (£=2.48-2.65), al- 
though in the latter case the underlying continuum was a 
factor of 1.5-4 brighter. The rather high value of the reflec- 
tion parameter (R~2) derived from the XMM-Newton spec- 
tra may be explained if the continuum source has fairly re- 
cently and/or temporarily entered a low-flux period, whereas 
the Compton reflection has remained high due to the time 
lag inherent in a component originating in material rela- 
tively distant from the central source of ionizing radiation. 
For example, the bulk of the Compton reflection may be pro- 
duced from the far side of a putative molecular torus in NGC 
4151, in which case the response time would be measured 
in years. The measured flux of the Compton reflection signal 
is a factor of 1.5iJ'g times higher than observed during the 
January 1999 BeppoSAX observation; this encompasses the 
possibility of little or no change in the Compton reflection 
signal over a two year interval. 



5.1 The profile of the iron Ka line 

Above we model the iron Ka emission observed in the EPIC 
spectra as a single, "narrow" (see §6.3) Gaussian component. 
The line energy is consistent with the fluorescence of neu- 
tral or near-neutral iron (taking into account up to 5 eV 
systematic error in the energy calibration). The line energy 
measured previously by ASCA was not sufficiently well con- 
strained to rule out the possibility that the iron Ka line is 
seen in purely scattered light (the first-order peak due to 
Compton scattering is lower by ~80 eV than the zeroth- 
order peak at 6.4 keV; Jourdain & Roques 1995). The iron 
Ka line energy measured in the XMM-Newton spectra rules 
out this possibility; the dominant line component is seen 
directly rather than as a result of a strong scattering of 
the iron Ka line in the extended medium above and below 
the accretion disk. Including a second Gaussian line com- 
ponent with a fixed energy of 6.4 keV and a significantly 
broader profile (e.g. a ~ 0.2 keV; Yaqoob et al. 2001) does 
not improve the overall fit; the upper limit on the flux in 
such a component is 1.0xl0~ 5 photons cm~ 2 s _1 . Fixing 
the energy of the broad Gaussian line at 6.2 keV (as found 
in previous studies - Yaqoob et al. 1995) results in a very 
weak line (the best-fit normalisation is an upper limit of 
5.1 xlO -6 photons cm -2 s _1 ) and does not improve the fit 
significantly (Ax 2 =2 for 2 d.o.f). Substituting the line emis- 
sion expected from an accretion disk (the disk-line model in 
XSPEC, Fabian et al. 1989), for the broad Gaussian compo- 
nent in this model, also results in a very weak line (best-fit 
normalisation is an upper limit of 2.75xl0 -6 photons cm -2 
s _1 ) that does not provide any improvement in the quality of 
the fit (Ax 2 =0.2 for 2 d.o.f). Replacing the two line profiles 
in these fits with a single line associated with emission from 
the accretion disk results in a similar y 2 to that obtained 
with a narrow line. However, in this case the derived inner 
radius of the accretion disk is r^, ~1000 r 9 (q~2, 8 «45°), 
which results in a disk-line profile that closely resembles a 
narrow Gaussian line profile. 

In summary, the high significance XMM-Newton spec- 
tra reveal only a narrow iron Kq line with no evidence for 
the presence of a strong relativistically broadened iron Kq 
feature. This is contrary to the earlier conclusions of Yaqoob 
et al. (1995) and Wang et al. (1999, 2001), who on the basis 
of ASCA spectra identified broad iron Kq features in the 
NGC 4151 spectrum. However, we note that the presence of 
a broad line was not required in the spectral modelling of 
the final ASCA long-look observation of NGC 4151 (Schurch 
& Warwick 2002; Takahashi et al. 2002). It remains possible 
that a relativistically broadened line appears in the X-ray 
spectrum of NGC 4151 only intermittently or for particular 
'states' of the source. Further monitoring of NGC 4151 with 
XMM-Newton and Chandra should help clarify whether this 
is the case. 



5.2 The iron Kq line flux 

The measured flux in the narrow iron Ka line is ~1.3xl0 -4 
photons cm -2 s~\ significantly lower than that measured 
in earlier ASCA, BeppoSAX and Ginga observations, where 
the line flux was more typically ~2.2xl0 -4 photons cm -2 
s _1 ; (Yaqoob & Warwick 1991; Yaqoob et al. 1993; Schurch 
& Warwick 2002). The present measurement is also some- 
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what lower than that reported in recent Chandra observa- 
tions (1.8±0.2xlCT 4 photons cm~ 2 s _1 ; Ogle ct al. 2000). 
Ogle et al. (2000) state that Chandra spatially resolves 
65±9% of the iron line emission, placing its origin in the 
ENLR at a distances of up to ~500 pc from the continuum 
source. Thus an iron Ka line flux of up to 1.2 x 10 -4 photons 
cm~ 2 s _1 could be attributable to the ENLR, a value com- 
parable to that measured in the current XMM-Newton ob- 
servations. The implication is that a component of the iron 
Ka line produced close to the central engine in NGC 4151 
almost completely disappeared over the ~10 month time 
interval between the Chandra and XMM-Newton observa- 
tions. However, we caution that this intriguing scenario may 
not be valid. The line intensity given by Ogle et al. (2000) 
combined with the exposure time and HETG effective area 
indicate ~160 photons in the line. The central ~l" of the 
cross- dispersion profile integrated over the energy interval 
containing the line contains the majority of the photons from 
both line and continuum (whether extended or not). To de- 
termine whether any of the line emission is extended to a 
given level of confidence, one needs to show that the cross- 
dispersion profile outside of the central ~l" lies above the 
background noise. Ogle et al. (2000) do not give a confi- 
dence level to their claim that the 1-3" interval either side 
of the peak of the cross-dispersion profile contains extended 
iron line emission, however a brief calculation of the cross- 
dispersion profile for the number of photons detected in the 
Chandra grating observation indicates that strongly signifi- 
cant (i.e. >4<J significance) spatially extended iron Ka line 
emission cannot be detected in this relatively short Chan- 
dra observation. If all the iron Ka line flux instead origi- 
nates close to the central engine, then the line flux from this 
vicinity would only need to decrease by ^35% between the 
XMM-Newton and Chandra observations. The intensity of 
the neutral iron Ka line in the XMM-Newton spectra is con- 
sistent with all the line flux originating in a nearly Compton 
thick molecular torus (the EQW with respect to the Compton 
reflection continuum, is 1.06 keV). The decrease in line flux 
observed between the BeppoSAX and XMM-Newton obser- 
vations is, however, in marked contrast with the behavior of 
the Compton reflection component, which remains at a sim- 
ilar flux-level or, if anything, is slightly stronger during the 
XMM-Newton observation than it was in the January 1999 
BeppoSAX observation. This decoupling of the iron Ka and 
the reflection continuum suggests that the additional iron Ka 
line flux observed in the earlier observations is formed in 
optically-thin, cool clouds located within about a light year 
of the central source (i.e. closer than the material respon- 
sible for the Compton reflection continuum). In this setting 
the decline of the iron Ka line flux presumably coincides with 
the central source entering a fairly prolonged low-flux period, 
sampled only latterly by the XMM-Newton observations. 

The variability constraints suggest that a significantly 
component of the iron Ka line may be associated with the 
complex absorbing medium, including the warm absorber. 
Schurch & Warwick (2002) note that the predicted range of 
iron ionization states in the warm medium, determined from 
a long-look ASCA observation, is Fe XVI-Fe XXI (for which 
the associated Ka line energy is 6.4-6.6 keV). Since the 
line energy measured by ASCA (Eku=6.395±0.012 keV), 
excludes much of this range, it was concluded that the bulk 
of the line flux derives from less strongly photoionized mate- 



rial located further from the nucleus than the warm medium. 
Comparison of the ASCA results with the present XMM- 
Newton observations suggests the narrow iron Ka line may 
typically be composed of a relatively constant component 
amounting to ^1.2xl0~ photons cm~ s~ plus a more 
variable component of comparable magnitude. In such a sit- 
uation the XMM-Newton line measurements clearly relate 
predominantly to the first of these two origins. The line 
equivalent width (^175 eV referenced to the observed power- 
law but ^60 eV for a more typical continuum level) and en- 
ergy derived from the XMM-Newton observations are con- 
sistent with an origin for the "constant" line component in 
reflection from a neutral of lightly photoionized medium lo- 
cated outside the immediate proximity of the nucleus i. e. r> 1 
pc. On the other hand the variable line emission component 
may still be associated with the warm absorbing medium pro- 
vided the dominant ionization state is closer to Fe XVI than 
Fe XXI (which could apply if our line of sight samples the 
high end of the ionization distribution in the warm medium). 
The typical equivalent width of ^50 eV for the variable com- 
ponent coupled with the estimate, NH,warm~2x 10 23 cmT 2 , 
then implies that the warm medium subtends an angle of 
^n/AFe steradians as viewed from the central source, where 
Afs is the iron abundance relative to solar. 

5.3 The intrinsic width of the narrow iron Ka line 

The global best-fit to the XMM-Newton data is obtained 
for an intrinsic line width of the "narrow" iron Ka line, 
o~Ka=33t_\ eV. However, continuing uncertainties in the 
CTI corrections for the current EPIC response matrices will 
impact on this measurement. An analysis restricted to the 
EPIC MOS spectra provides a good check on the value of 
the line width, despite the marginally worse spectral reso- 
lution at 6.4 keV, since the calibration is relatively better 
known for the MOS cameras compared with the PN camera. 
This analysis used epoch-specific response matrices, incor- 
porating time dependent CTI corrections (S. Sembay, pri- 
vate communication) and yielded an intrinsic line width of 
o~Ka=51~)L 7 7 eV (Note: instrument response matrices with a 
non-standard 1 eV energy scale were used for this analysis). 
By way of comparison the effective MOS instrument reso- 
lution at the observation epoch was a~70 eV (FWHM~165 
eV). At face value the EPIC spectra suggest an intrinsic 
width of between 3000-6000 km s' 1 (FWHM) for the iron 
Ka line. In contrast, Ogle et al. (2000) report, on the ba- 
sis of Chandra HETG observations, that "the narrow core 
of this line is unresolved, with FWHM=1800±200 km s' 1 ". 
Clearly this cautions against over-interpreting the EPIC line 
widths, particularly since the measurements are close to the 
limit of the present calibration. 

What might be the origin of such line broadening if con- 
firmed? One consideration is that the iron Ka line is in fact 
a doublet with components at 6.404 (Ka\) and 6.391 (Ka?,) 
keV. However, we find that modelling the iron Ka emission 
feature with two Gaussian components with a fixed 13 eV 
energy separation and flux ratio of 2:1 (as per the branch- 
ing ratio) has little impact on the above results. A further 
possibility is that the Ka line appears broadened because it 
represents a blend of ionization states. However, the mea- 
sured line energy is fully consistent with that of neutral iron 
at 6.400 keV (representing the centroid of the Ka\ and Ka? 



© 0000 RAS, MNRAS 000, 000-000 



6 N.J. Schurch, R.S. Warwick, R.E. Griffiths S. S embay 



components) and we find that adding non-neutral line com- 
ponents fails to replicate the low energy wing of the line. 
Finally, the line broadening could be due to Doppler shifts, 
which would place the origin of the Ka line observed by 
XMM-Newton squarely in the broad line region of the galaxy. 
Unfortunately this description is then at odds with our ear- 
lier conclusion that the relatively "constant" component of 
the line emission (observed by XMM-Newton) probably orig- 
inates in reflection from lightly photoionized matter beyond 
pc from the nucleus. 

5.4 The additional iron edge - a clue to the iron 
abundance. 

The best-fit spectral model includes iron edge features im- 
printed on the direct continuum by the neutral and warm 
absorbers and also arising in the Compton reflection com- 
ponent. In each case solar metal abundances are assumed. 
However, as detailed earlier the high signal-to-noise EPIC 
spectra also require the inclusion of an additional absorption 
above 7 keV suggestive of an overabundance of iron in the 
absorbing and/or reflecting media. Modelling the hard X-ray 
spectrum in the region of the iron edge with a simple model 
incorporating a power-law continuum and a single absorp- 
tion edge yields an edge energy of 7.11_q 01 keV and a total 
optical depth of TK,Edge~0.23^Q'g2- The depth of this edge 
overpredicts the observed iron line flux by a factor of ^2.5, 
implying a covering fraction for the absorbing material of 
^4-0%. The derived edge energy is consistent with absorption 
by neutral iron (with the upper limit of 7.13 keV, including 
a systematic uncertainty of 5 eV, formally excluding even 
Fe II). This rules out an association of the excess absorp- 
tion with the warm absorber, for which the predicted range 
of edge energy is 7.8-8.3 keV (Fe XVI-Fe XXI). If we inter- 
pret the observed optical depth in the edge (TK,Ed g e~0.12) as 
due to an extra-solar abundance of iron in the cold absorber 
we obtain, A Fe ^3.6, whereas its association with the Comp- 
ton reflector would imply Af^2.3. Alternatively if we allow 
an overabundance in both the cold absorber and the reflector 
then Af^2 is required. Previous studies have also reported 
evidence for an overabundance of iron in NGC 4151, for ex- 
ample Yaqoob et al. 1993 quote a canonical value of Af^2.5 
times solar, based on the average of many Ginga measure- 
ments. 



6 THE CONCLUSIONS 

Recent observations with the EPIC cameras on XMM- 
Newton have provided a very high signal-to-noise measure- 
ment of the hard (2.5-12 keV) X-ray spectrum of the archety- 
pal Seyfert galaxy NGC 4151. Remarkably, we find that the 
spectral "template" model developed by Schurch & Warwick 
(2002) to fit earlier BeppoS AX and ASCA observations, also 
provides an excellent description of the XMM-Newton data, 
with only minor modifications. 

This supports the view that the complex spectral vari- 
ability exhibited by NGC 4151 over timescales longer than a 
few days is the product of changes in the source, driven pre- 
dominantly by the level of the underlying continuum. The 
ingredients of the spectral model are surprisingly simple. An 
underlying power-law continuum is absorbed below ~5 keV 



by a combination of warm and cold gas along the line-oj '-sight 
but supplemented at higher energy by Compton-reflection. 
The parameters specifying the continuum slope, the column 
densities and the Compton reflection flux remain essentially 
constant. There is no requirement for a relativistically broad- 
ened iron Ka line feature, but the intensity of the nar- 
row component definitely varies on long (years) timescales. 
There is also a strong suggestion that either the cold absorb- 
ing and/or reflecting media have an iron abundance that is 
twice solar. 

Further development of this picture of the hard X- 
ray properties of NGC 4151 will be possible through future 
approved XMM-Newton observations aimed at monitoring 
spectral changes in the 2.5-12 keV band. With the addition 
of a simultaneous INTEGRAL observation it will be possible 
to test the assumption that it is the level of the broadband 
hard X-ray continuum which governs such changes. 
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